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Long range interactions between like homonuclear alkali metal diatoms 

Jason N. Byrd/'E3 Robin Cote/ and John A. Montgomery, Jr.^ 

Department of Physics, University of Connecticut, Storrs, CT 06269 

Long range electrostatic and van der Waals coefficients up to terms of order R^^ have been evaluated by the 
sum over states method using ah initio and time dependent density functional theory. We employ several 
widely used density functionals and systematically investigate the convergence of the calculated results with 
basis set size. Static electric moments and polarizabilities up to octopole order are also calculated. We present 
values for IA2 through K2 which are in good agreement with existing values, in addition to new results for 
Rb2 and Cs2. Interaction potential curves calculated from these results are shown to agree well with high 
level ab initio theory. 



I. INTRODUCTION 

The long range interactions between atom-diatom and 
diatom-diatom systems at low densities and cold temper- 
atures is important in many areas of atomic and molec- 
ular physics. Applications may be found in precision 
spectroscopyji condensed matter physics^ as well as in 
the studyi"— and control^ of chemical reactions. In the 
study of low density diatomic collisions, interactions at 
long range play a critical role in the reaction process ii 
The weakness of the long range intermolecular interac- 
tion compared to the chemical bond and the range of 
nuclear coordinates necessary to model a long range po- 
tential energy surface (PES) suggests that it is advanta- 
geous to consider modeling the intermolecular potential 
in other ways than brute force ah initio quantum chemi- 
cal calculations. 

In the limit where the wavefunction overlap between 
two interacting monomers is negligible, the interaction 
potential between them can be expanded in a van der 
Waals series, where the interaction energy can be broken 
into three distinct components 



Eint — Eel + -E'ind 



E. 



disp ■ 



(1) 



Here Egi, Eind and -Edisp are the permanent electrostatic, 
induction (permanent-induced electrostatic) and disper- 
sion energies and can each be written in an asymptotic 
series 



Elr — ^ VnR 



(2) 



These contributions to the long range interaction energy 
can be calculated in several ways. In this work we expand 
the interaction operator in terms of multipole moments)^ 
and then calculate the interaction energy using first and 
second order Rayleigh-Schrodinger perturbation theory. 

Isotropic dispersion interactions between Na2 and K2 
pairs have been previously studied using the London 
formulaiS and time dependent density functional the- 
ory (TD-DFT)^iii Spelsberg et alM have calculated 



van der Waals coefficients including anisotropic correc- 
tions using valence full configuration interaction (VFCI) 
theory for Li2, Na2 and K2. Also recent work by 
Kotochigovai^ has calculated the isotropic van der Waals 
interaction for KRb with leading order anisotropic cor- 
rections. In this paper we present results for the van 
der Waals interactions between pairs of ground state 
homonuclear alkali metal diatoms, through order R~^, 
including anisotropic corrections. In Sec. |TI]an outline of 
the sum over states method of calculating van der Waals 
coefficients is presented, followed by a discussion of the 
electronic structure calculations done to obtain the van 
der Waals coefficients in Sec. IIIIl We conclude in Sec. 
II VI with a discussion of our results. The extension to het- 
eronuclear alkali diatomic van der Waals interactions will 
be presented in subsequent publications. 



II. ANISOTROPIC LONG-RANGE 
INTERACTIONS 

In the Born-Oppcnheimcr approximation, the interac- 
tion between two linear molecules can be expanded in a 



FIG. 1. Basis set convergence of the octopole static polar- 
izability values for Na2 evaluated using several ab initio and 
DFT methods. 
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Dispersion 
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8.304(3) 


4.679(3) 


4.094(3) 


3.697(3) 


3.654(3) 


4.216(3) 


4.073(3) 
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3.091(2) 
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5.072(2) 
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5.862(1) 
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1.401(5) 


1.086(5) 


1.073(5) 


4.984(4) 


6.712(4) 


7.579(4) 


7.459(4) 
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9.006(4) 


6.887(4) 


6.834(4) 


3.237(4) 


4.335(4) 


4.909(4) 


4.802(4) 
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2.511(4) 


1.390(4) 


1.784(4) 


2.067(4) 
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8 2 2 1 


-7.153(3) 


-4.849(3) 


-5.022(3) 


-2.779(3) 


-3.569(3) 


-4.134(3) 


-3.855(3) 
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1.788(3) 
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6.003(4) 
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2.355(4) 


2.728(4) 


2.545(4) 


8 4 2 1 


-6.438(3) 


-4.364(3) 


-4.520(3) 


-2.501(3) 


-3.212(3) 


-3.720(3) 


-3.470(3) 


8 4 2 2 


3.576(2) 


2.425(2) 


2.511(2) 


1.390(2) 


1.784(2) 


2.067(2) 


1.928(2) 


%RMS 


22.1 


7.3 


8.0 


5.5 


8.6 


3.0 


4.3 



(2) 

TABLE I. Unique calculated dispersion and induction van der Waals coefficients, W^j^_^];^^i^.i, for ground state Na2 using selected 

ah initio and DFT methods. The RMS deviations are relative to previous theoretical results.— All calculations are performed 
at the experimental equilibrium bond length. 



complete angular basis as 

£'int(ri,f2,R) = ^ EL^L2L{R)AL^L2L{fi,f2,R), 

L1L2L 

(3) 

where fi = {6i,(j)i) are the molecular orientations and 
R = {R, 6, (j)) defines the vector between the molecu- 
lar centers. As shown by Mulder, van der Avoird, and 
Wormer,^ if we choose coordinates so that R is oriented 
along the z axis, the angular functions may be written as 

min{Li L2) 

AL,L2Lifl,f2,R) = Vl[l2L 
M=0 

X (cos ei)Pll{cOS 62) COS Af(0i - 02), (4) 

where 

vI[l2L = (-1)*'(2 - 5Mfl){LiM- L2 - Mm 

■ {L,~M)\{L2-M)\ y'^ 
_{Li+M)\{L2 + M)\\ ' ^ ' 



{LiM;L2 — M\LO) is a Clebsch- Gordon coefficient, and 
P]^{cos9) is an associated Legendre polynomial. Be- 
cause the interaction energy is rotationally invariant, 
it may be expanded in terms of multipole operators 
Qim = J2i^i'^i^im{ri), where Cimifi) is a Racah spher- 
ical harmonic.i^ Using first- and second-order perturba- 
tion theory, Mulder et ali^^ express the coefficients 
EliL2l{R) in terms of the separated monomer transi- 
tion moments. When this is done, the interaction energy 
may be written as 

min{Li ,L2) 

E,M0i,ci>ue2,^2)^Y. E ^lIlIAR) 

Pll (cos ei)Pll{cOS 62) COS Af(0i - 02). (6) 

Here the first-order contribution 

^LiLmC-^) = ^nLiL2MR^''^Li+L2 + l,n (7) 
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785.8 


448.3 
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7.956 


789.7 
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CS2 


B3PW91 


8.78 


1096.0 


577.8 




PBEO 


8.78 


1032.9 


568.3 




Ref. 18 


8.78 


1012.2 


509.0 



TABLE II. Static dipole polarizability values for the X^Ti^ 
homonuclear alkali diatoms calculated at the experimental 
equilibrium bond length. 



(1) 



nLiL2M 



{Li+L2)l 



[Li + M)\{L2 + M)\ 
X (0i|Ql,o|0i)(02|Ql,o|02) (8) 



is due to the electrostatic interaction and the second- 
order terms 



V, 



(2) 



L1L2M 



(i?) 



^ 0h+l[+l2+l'2+2,n 



(9) 



a 



kik2 



(10) 

contain contributions from dispersion and induction. The 
implies that A:i + fe2 ^ and Ck^ is the energy of the 

/ci'th state. The ClilHi coefficient is a scalar coupling 
term given hy^ 

ClS = (-l)'^+'H(2ii + 1)!(2L2 + 1)!)^/' 

■(2;i + 2;2 + i)!(2;i + 2/^ + i)!T'/' 



L2L 



x(Zi+Z20;;l + ZiO|LO)<( h 1'2 L2) , (11) 

h + h l'i+l'2 L 



B3PW91 

PBEO 

SAOP 

London_^ 

Valence~FCI'^ 



Li2 

2.649(3) 
2.593(3) 



2.730(3) 



Naa 

4.216(3) 

4.073(3) 

4.460(31^ 

4.374(3) 

4.181(3) 



K2 

1.242(4) 

1.185(4) 

1.106(4)^ 

1.186(4) 

1.039(4) 



Ref. [TO 
Ref. [O 
" Ref. la 
Ref. [12 

TABLE III. Comparison of theoretical isotropic Va van der 
Waals coefficients. 



defined as 

T°jf£^ =^(0,|gi,™|fc,)(A:,|Q,,_„jO,)(;,TO;?^ -m|i,0) 

m 

(12) 

where the indices ki go over ground and excited states. 
These coupled transition moments also transform the 
same way as the electrostatic moments {Qlm)o- It 
should be noted that from Eq. [TUl contributions such 
as 



k2^0 



2^02^2 
^2^o^2 



eo2 



+ (1^2) 



(13) 



are associated with (pi)^a2 + (1 2) type induction 
terms. Contributions to Eq. [TOlsuch as 



rpOiki rp02k2 



fe2#0 



£02 



(14) 



are associated with aiQ;2 type dispersion terms and can 
can be related to a Casimir-Polder integral over complex 
frequencies of coupled dynamic polarizabilities (see Ap- 
pendix 



It is convenient in practice to collect all terms with the 
same i?-dependence in Eq. [9] into a single expression 



V, 



(2) 



L1L2M 



{R) = - 



(2) 



nL\L2M 



(15) 



The dispersion and induction contributions to Eq. I15l are 
calculated separately as 

^(2) _ ^r(2,DIS) ^(2JND) , . 



The coupled transition moment for each monomer is and presented in tables |T] and |Vl 



3 



n Li L2 M 


Li2 


Na2 


K2 


Rb2 


CS2 


Q20 


10.74 


10.52 


15.68 


16.06 


27.85 


5220 


6.921(2) 


6.640(2) 


1.475(3) 


1.548(3) 


1.912(2) 


5221 


-1.538(2) 


-1.476(2) 


-3.278(2) 


-3.439(2) 


-4.248(2) 


5222 


9.612(0) 


9.223(0) 


2.049(1) 


2.149(1) 


2.655(1) 



TABLE IV. Electrostatic van der Waals coefficients, W, 



(1) 



for the ground state alkali dimers Li2, Na2, K2, Rb2 and Cs2 



calculated at the CCSD(T) level of theory using the finite field method. 



III. ELECTRONIC STRUCTURE 
CALCULATIONS 

Transition moment calculations for use in the eval- 
uation of Eq. [T^l were done for the homonuclear al- 
kali metal diatoms Li2, Na2, K2, Rb2, and Cs2 in the 
X^E+ ground state using a locally modified version of 
the GAMESSi-!- - quantum chemistry program package. 
All calculations have been performed at the experimen- 
tal equilibrium bond distances. The number of excited 
states included in the sum in Eq. [10] is taken to be the 
total number of virtual orbitals present. Tests were per- 
formed to ensure convergence of the calculated TD-DFT 
transition moments with respect to the grid size. The 
grid employed in production calculations uses 155 radial 
points for all atoms, and prunes from a Lebedev grid 
whose largest size is 974, thus using about 71,000 grid 
points/atom (the JANS=2 grid in GAMESS). 

To provide consistent results for all the alkali metals, 
the Karlsruhe def2-QZVPP— basis sets were chosen for 
this work. The def2 basis sets are available for almost 
the entire Periodic Table and are well-known for their 
robustness and their excellent cost-to-performance ratio 
in large-scale Hartree-Fock and density functional theory 
(DPT) calculations. The def2-QZVPP basis sets for the 
alkali metals contain spdf basis functions with two po- 
larization functions. For the row five and lower atoms, 
the inner-core electrons are replaced by an effective core 
potential (ECP) to reduce the number of electrons in- 
cluded in the correlation treatment and to account for 
scalar-relativistic effects. The ECP-28^ and ECP-462^ 
Stuttgart pseudopotentials were used for Rb and Cs re- 
spectively. These ECP definitions leave the sub-valence 
s and p electrons free, which are known to contribute the 
most to core- valence correlation energy^'^ in alkaline sys- 
tems. In the transition moment calculations presented 
here, core-valence correlation is accounted for implicitly 
within the DFT formalism. 

When examining the effects of basis set convergence of 
the van der Waals coefficients it is convenient to first look 
at the convergence of the multipole polarizability defined 
as 



2^ (0|g;™|fc)(fc|g;'-™|0} 



(17) 



The octopole polarizability terms 01300 are the most sensi- 



tive to the effects of higher angular momentum functions 
in the basis set, and were therefore used to asses the 
convergence of our results with respect to basis set size. 
Test calculations including a sequence of g basis functions 
ranging from valence to diffuse indicated that there was 
no need to add higher angular momentum functions to 
the def2 basis sets. The effects of adding additional sets 
of diffuse even tempered spdf basis functions can be seen 
in Fig. [TJ It was found that two additional sets of diffuse 
functions was necessary to achieve accurate results for 
our test case of Na2. For a subset of the DFT function- 
al considered here, the addition of a third set of diffuse 
functions was investigated and found to contribute little 
to the calculated polarizabilities. The def2-QZVPP basis 
sets for K, Rb and Gs include an additional spdf diffuse 
function by definition and so only required a single ad- 
ditional set of diffuse basis functions. We note that the 
dipole polarizability, and hence the Ve van der Waals co- 
efficients, are essentially converged using the unmodified 
def2 basis set. A similar procedure was recently used by 
Rappoport and Furche^^ to optimize the Karlsruhe def2 
basis sets for molecular property calculations. 

Using the double augmented basis set (d-aug-def2- 
QZVPP) described above two ab initio methods were 
considered for our test case Na2, configuration interaction 
singles^^''^" (CIS) and time dependent Hartree-Fock^Si^ 
(TD-HF) theory. Additionally we used time dependent 
density functional theory (TD-DFT)^i^ to evaluate ex- 
citation energies and multipole transition moments. The 
results of the ab initio and TD-DFT calculations are com- 
pared to those of Spelsberg and Meyer— in Table |T1 De- 
spite the success^ in describing the van der Waals inter- 
action of II2 , CIS significantly over estimates the disper- 
sion and induction van der Waals interaction. This is not 
unexpected as CIS is known not to obey the Thomas- 
Reiche-Kuhn dipole sum rule^^"— because it overesti- 
mates the transition moments. An additional source of 
error in CIS is the overestimation of excitation energies 
that results from the use of canonical Hartree-Fock vir- 
tual orbital energies. The overall performance of TD- 
DFT compared to the results of Spelsberg and Meyer- 
appears quite good. Upon examining the results of Table 
|T]we have chosen the B3PW91 and PBEO functionals for 
use in the rest of this work due to the consistent accuracy 
for both dispersion and induction interactions. 

The electric quadrupole moment used to calculate 
V5 was evaluated using the finite field method using 



4 



CCSD(T). Core- valence contributions are accounted for 
by including the inner valence s and p electrons in the cor- 
relation treatment. For consistency the n-def2-QZVPP 
basis sets defined above were used. Comparisons with 
the results of Harrison and Lawson'^'^ for Li2 , Na2 and K2 
show that this choice of basis set provides reliable results. 
In Table HV] we show our calculated quadrupole moment 
and Vj; van der Waals coefficients for the homonuclear 
alkali diatoms including new results for Rb2 and Cs2. 



IV. NUMERICAL RESULTS AND 
DISCUSSION 

As noted above, the Karlsruhe def2 basis sets replace 
the inner shell electrons of rubidium and cesium with an 
effective core potential. In order to assess the effect of re- 
placing the inner-core electrons in heavy atoms by a pseu- 
dopotential, we have performed calculations of the K2 
van der Waals coefficients using the Stuttgart ECP and 
basis set and made comparisons with our all-electron 
calculations using the def2 basis set. The Stuttgart sp 
ECP basis set was uncontracted, and to ensure that the 
basis set was fully saturated the most diffuse five p ex- 
ponents were added as d functions and the subsequent 
most diffuse four d exponents were added as / functions. 
Using this basis and ECP, it was found that the differ- 
ence between the ECP and all-electron def2 results were 
negligible. 

An important consideration for any van der Waals 
study is that of the dipole polarizability of the individual 
molecules as the leading dispersion interaction term Vg 
is proportional to the product of polarizabilities. Using 
our n-aug-dcf2-QZVPP basis sets with the B3PW91 and 
PBEO DFT functionals we have calculated the dipole, 
quadrupole and octopole polarizabilities for all the alkali 
diatoms considered in this work. In Table HI] the dipole 
polarizabilities are presented with several other existing 
theoretical results in comparison. The calculated higher 
order polarizabilities can be found in the supplemental 
materialj56 Because the polarizability is very sensitive to 
the diatomic bond distance near equilibriumjiS we have 
chosen to perform all calculations at the experimental 
equilibrium bond length, rather than at the optimized 
geometry for each different level of theory. As a con- 
sequence, there will be small deviations in the parallel 
components of the calculated polarizability at a particu- 
lar level of theory compared to what would be found at 
the optimized geometry. Therefore, our calculated values 
for the perpendicular dipole polarizability (am) agrees 
better with the existing theoretical results than the par- 
allel polarizability (ano). For the higher order polariz- 
abilities, the only available values come from Spelsberg 
et alJ^. The total polarizability RMS error for Li2, Na2 
and K2 is 8.7%, 6.7% and 9.5% and 7.5%, 5.1% and 6.9% 
for B3PW91 and PBEO respectively. 

Often neglected, the leading order term in the van 
der Waals expansion for homonuclear alkali diatoms is 



1 , ' 1 ' 1 


' 1 ' 1 ' _ 




0-0 CCSD(T)-F12a/QZVPP 






-- -c;"/r' 






— Cj/r' - Cj'^/r' 






— vdW TD-DFT PBEO 








' \ ! . 1 , 1 


1,1,' 



R (a.u.) 

FIG. 2. Comparison of the van der Waals surface of co-linear 
Na2 with both the London isotropic Ve approximation^ and 
a fully ab initio curve. 



the electrostatic V5 quadrupole-quadrupole interaction. 
This leads to a repulsive interaction potential for a sig- 
nificant portion of the interaction phase space with a 
long range barrier; for the case of Na2 this barrier has 
a maximum height of ^ 13cm~^ for the co-linear orien- 
tation [Oi —02 — 0). There are however orientations 
of approaching pairs of molecules for which no barri- 
ers are found. For example, this was seen for Na2 and 
K2 by Zemke et aM- as well as for the related systems 
of K2-I-Rb2^ and Rb2-I-Cs2^ where the lowest energy 
reaction path was calculated ab initio and indeed was 
found to be barrier less. To illustrate the effect of in- 
cluding the V5 electrostatic term as well as the higher 
order Vg terms, the results of the London approximation 
given by Zcmkc et al.^ are compared in Fig. [5] to the 
TD-DFT/PBEO induction and dispersion coefficients re- 
ported above for the co-linear orientation. As a compar- 
ison, the potential curve was also calculated using the 
CCSD(T)-F12a (explicitly correlated CCSD(T)) level of 
theory^i^. As can be seen in Fig. [31 the inclusion of 
V5 greatly improves the long range representation of the 
potential, while the higher order van der Waals terms are 
necessary to describe the interaction as the molecule sep- 
aration decreases. With the inclusion of the higher order 
terms reported in this work, the van der Waals surface 
matches the barrier height and molecule separation to a 
few percent. 

In a cold ( J = 0) ensemble with no external field the 
rotationally averaged electrostatic V5 van der Waals co- 
efficient evaluates to zero. For this case the leading order 
term in the van der Waals series is the isotropic Vg dis- 
persion interaction. We compare our calculated Vg coeffi- 
cients to existing theoretical results in Table IIIII For the 
case of Li2 it can be seen that B3PW91 performs well, 
which is also true from Table U for the case of Na2 . For 
the heavier K2 we find that PBEO performs marginally 
better for Vg. In Table IVl all the calculated van der Waals 
coefficients including all anisotropic contributions for Li2 
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6.032(4) 


3.075(5) 


3.447(5) 


3.390(5) 


4.223(5) 


4.199(5) 


6.083(5) 


6.174(5) 


8 2 


3.837(4) 


1.236(5) 


2.254(5) 


2.196(5) 


2.779(5) 


2.736(5) 


4.036(5) 


4.057(5) 


8 2 2 


1.393(4) 


2.717(4) 


1.016(5) 


9.269(4) 


1.308(5) 


1.203(5) 


2.000(5) 


1.889(5) 


8 2 2 1 


-2.787(3) 


-3.771(3) 


-2.032(4) 


-1.854(4) 


-2.615(4) 


-2.407(4) 


-4.000(4) 


-3.778(4) 


8 2 2 2 


6.967(2) 


1.629(2) 


5.080(3) 


4.634(3) 


6.538(3) 


6.016(3) 


1.000(4) 


9.445(3) 


8 4 


2.585(4) 


1.698(3) 


1.477(5) 


1.453(5) 


1.810(5) 


1.800(5) 


2.607(5) 


2.646(5) 


8 4 2 


1.839(4) 


1.284(3) 


1.341(5) 


1.224(5) 


1.726(5) 


1.588(5) 


2.640(5) 


2.494(5) 


8 4 2 1 


-2.508(3) 


-1.680(2) 


-1.829(4) 


-1.668(4) 


-2.354(4) 


-2.166(4) 


-3.600(4) 


-3.400(4) 


8 4 2 2 


1.393(2) 


8.254(0) 


1.016(3) 


9.269(2) 


1.308(3) 


1.203(3) 


2.000(3) 


1.889(3) 


%RMS 


4.5 


5.4 


7.2 


5.9 











TABLE V. Unique van der Waals coefficients, W^^'^^^a/i ^'^i' the ground state alkali dimers Li2, K2 and new results for Rb2 and Cs2 calculated at the TD-DFT level 
of theory using the B3PW91 and PBEO functionals. The RMS deviations are relative to previous theoretical resultai^ where applicable. 



and K2 are presented. To compare with the results of 
Spelsberg et al^ it is necessary to recouple their inter- 
action coefficients from the Li,L2,L coupling scheme to 
the Li,L2,M coupling scheme used here by the following 

VL1L2A/ = ^Vl[l2lVliL2L- (18) 

L 

After recoupling we find our van der Waals coefficients 
to be in good agreement with those of Spelsberg et al. 
Also in Table |V] are our new results for Rb2 and Cs2 . 
Using the aug-def2-QZVPP basis set and ECP, our cal- 
culations for Rb2 and Cs2 contain comparable numbers 
of electrons to the all-electron calculations for Na2, and 
therefore we expect our van der Waals coefficients for 
these alkali diatoms to be comparably accurate. 



trap. Coupling with an external field can lead to a richer 
interaction phase space beyond that of the isotropic dis- 
persion interaction. 
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Appendix A: Relationship of the sum over states 
method to the Casimir-Polder formulation 



V. CONCLUSION 

In summary, we have used ab initio and TD-DFT 
methods to calculate excitation energies and transition 
moments for the homonuclear alkali diatoms Li2, Na2, 
K2, and Cs2. These energies and moments were used 
in a sum over states calculation of the electrostatic, in- 
duction and dispersion van der Waals coefficients as well 
as static polarizabilities up to third order. The Karl- 
sruhe def2-QZVPP basis sets were augmented with dif- 
fuse functions to ensure accurate calculation of the multi- 
pole polarizabilities and van der Waals coefficients. Cal- 
culated dipole polarizabilities were found to be consis- 
tent with other theoretical results, while noting that the 
polarizability is sensitive to the diatomic bond length.^® 
Van der Waals results for the lighter systems are in good 
agreement with the existing literature, and suggest errors 
of a few percent for the new results for Rb2 and Cs2. As 
such the long range interaction of each alkali diatomic 
pair can be characterized by a van der Waals series of 
the form 



The sum over states formulation of Eq. [TU] is not the 
only way to obtain dispersion interactions as described 
in Eq. 1141 We define the coupled dynamic polarizability 
as 



(Al) 



where Ae^; — ek^ — eo; is the fc^'th excitation energy. 
Using the identity 



1 



ab 



a + b TT Jo (a2 -Hw2)(62 +a;2) 
then Eq. [T3] can be rewritten as 



duj (A2) 



- / 4[Htu;)af^{tu;)doj (A3) 
Jo 

which is the well known Casimir-Polder integral for dis- 
persion interactions between two molecules. 



nLiL2M 

X (COS(01))F*^(COS(02)) COs(A/(</)i - </)2)) (19) 



where the 



(1.2) 

LjLnM 



coefficients can be taken directly from 



Tables U and |Vl We provide a sample FORTRAN pro- 
gram for evaluating Eq. [19] in the supplimental material 
for the case of Na2 as described in Fig. [21 The current 
study suggests that the calculation of the van der Waals 
expansion coefficients can provide an accurate model of 
the long range potential surface in alkali dimers. 

It is noted that the leading order van der Waals coeffi- 
cient is the electrostatic V5 quadrupole-quadrupole inter- 
action. This term provides a molecular frame repulsive 
interaction at long range for a portion of the diatom- 
diatom phase space. This repulsive interaction will lead 
to a stabilization of the molecules when confined in a 
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